﻿Modelare si simulare Seminar 1 SEMINAR 1 1 Experiment aleator Un fenomen a carui evolutie difera semnificativ atunci cand este repetat in aceleasi conditii se numeste experiment aleator Specificarea experimentului aleator consta in stabilirea procedurii de desfasurare a sa, precum si a setului de masuratori si observatii care-l insotesc Modelarea experimentului aleator este posibila daca evolutia sa se caracterizeaza prin regularitate statistica (expl mediile ) Rezultatul indivizibil obtinut in urma desfasurarii unui experimen t se numeste realizare sau esantion punct ( ζ ) Spatiul realizarilor sau spatiul de esantionare (S) = multimea tuturor realizarilor posibile I Spatiu de esantionare discret = multimea realizarilor unui experiment care e finita sau infint numarabila II Spatiu de esantionare continuu = multimea realizarilor experimentului e infinit nenumarabila Se numeste eveniment (A) o submultime a spatiului de esantionare ¾ daca AS= eveniment sigur ⇒ ¾ daca A=φ ⇒ eveniment imposibil sau eveniment nul ¾ in spatiul de esantionare discret un eveniment ce consta dintr-o singura realizare = eveniment elementar Fiind multimi, evenimentele se pot combina cu ajutorul operatiilor cu multimi (, ,)∪∩ sau se pot compara cu ajutorul relatiilor dintre multimi (incluziune, egalitate) Aplicatia 1 Fie urmatoarele experimente: E1: la receptie se detecteaza trei simboluri; se noteaza secventa binara corespunzatoare; E2: la receptie se detecteaza trei simboluri; se numara simbolurile de 1 (receptionate) continute; 1Modelare si simulare Seminar 1 OBS Evenimentele E1 si E2 se aseamana prin procedura dar difera prin masuratori E3: un nod de comutatie deserveste N terminale; se noteaza nr de terminale active la un moment dat; E4: un bloc de informatie este repetat la transmisie pana cand primeste un mesaj de confirmare; se numara repetarile; E5: la un centru de distribuire sosesc mesaje; se masoara durata intre doua sosiri succesive; E6: o componenta e introdusa intr-un montaj ce se pune in functiune; se masoara timpul de functionare al componentei; E7: la receptie se urmareste semnalul din linie; se masoara si se noteaza tensiunea la momentul t1; E8: la receptie se urmareste semnalul din linie; se masoara si se noteaza tensiunile la momentele t1 si t2; i) Sa se determine spatiile de esantionare corespunzatoare; ii) Considerand experimentul E7: se definesc urmatoarele evenimente: A - “ amplitudinea tensiunii e mai mare de 10V” B – “tensiunea e mai mica de -5V” C – “tensiunea e pozitiva” Sa se exprime in limbaj matematic evenimentele: CC A, B, C, AB∪, AB∩, C()ABC∪∩, , ABC∩∩, ()AB∪ C unde ()reprezinta complementul evenimentului Rezolvare: i) 000, 001, 010, 100, 011, 110, 101, 111S = {} 1 0, 1, 2, 3S = {} 2 { 0, 1, ,SN=… } 3 { 1, 2, 3,S =… } 4 0[0,SS tt== ≥=∞) { } 56 {(,Su u=−∞ =−∞−∪+∞ {} { 5(,5Buu= =+)∞ { } {; 5 10 ( , 5) (10, )AB uu sauu∪ = = −∞ − ∪ +∞ } 2Modelare si simulare Seminar 1 {;10(,1AB uu∩= = +∞ () {;10}(10,ABCuu∪∩) ()ABC∩∩=φ C {}( ) ; 5 10 [ 5,10]AB u u∪=−≤≤=− Frecventa relativa (f) de aparitie a unui eveniment A pe parcursul repetarii de n A ori a unui experiment este data de relatia: ()Nn A f= A n ()Nn- frecventa absoluta a evenimentului A A OBS: In conditii de regularitate statistica raportul pe termen lung (pt un numar foarte mare de repetari): ()fnp= limAA(1) n→∞ probabilitate de aparitie a evenimentului A Legatura dintre aceasta marime abstracta probabilitatea si marimea intuitiva, frecventa relativa este oferita de doua teoreme importante: 1 Legea slaba a numerelor mari ()Nn A Exprima probabilitatea ca raportul p→ sa tinda catre 1 cand n→∞ A n 2 Teorema limita centrala “Masoara abaterea dintre modelul ales (in cazul de fata valoarea p) si frecventa A relativa ()fn obtinuta in urma unui numar finit de incercari A Utilizarea relatiei (1) ca definitie a probabilitatii genereaza dificultati in dezvoltarea unei teorii a probabilitatilor, deoarece: nu este clar in ce sens exista limita ; un experiment nu poate fi repetat la infinit ; prezinta inconsistenta in cazul unui eveniment cu aparitie rara De aceea are sens dezvoltarea unei teorii abstracte a probabilităţilor, plecând de la premisele: 3Modelare si simulare Seminar 1 I Un experiment aleator a fost definit prin procedura de desfasurare, masuratorile si observatiile ce se efectueaza iar spatiul de esantionare S a fost identificat II S-a specificat (stabilit) o familie de evenimente Ω astfel incat: 1 S∈Ω si ; φ∈Ω C Ω→ ∈Ω 2 ()AA∀∈ 3 AA∈Ω → ∪ ∈Ω nNn∈ nN∈ (Ω este un tribut pe S, iar cuplul (S, Ω) se numeste spatiu masurabil) III Fiecarui eveniment A i s-a alocat (atasat) un numar (notat []PA) numit probabilitate astfel incat urmatoarele axiome sunt satisfacute: Axioma 1 0[]PA≤ Axioma 2 [] 1PS = Axioma 3 Fie A, B doua evenimente astfel incat: AB∩ =φ, atunci: [ ] [] []PAB PA PB∪= + ’ Axioma 3 Fie ,,AA… o secventa de evenimente astfel incat: 12 ∞ ∞ ,() ,AAij∩ =φ ∀ ≠atunci [] []PAP=A ∪ ijkk∑ 11kk= = ]([P… ca o masura a probabilitatii este o aplicatie a lui Ω in R, iar tripletul (, , )SRΩ se numeste spatiu probabilizat Un model matematic este probabilistic atunci cand conditiile de efectuare a experimentului vizat (ce prezinta regularitate statistica) determina exact probabilitatile de aparitie a realizarilor sale ¾ Construirea unui model probabilistic consta in particularizarea premizelor I, II, si III, la cazul, situatia reala studiata ¾ Ca orice modelare, calculul probabilitatilor este o tehnica matematica de descriere a fenomenelor aleatorii Validitatea sa intr-o aplicatie concreta nu poate fi infaptuita decat cu ajutorul metodelor statistice de interpretatre a observatiilor ¾ Regula prin care unui eveniment A al unui experiment E i se ataseaza o probabilitate []PA se numeste lege de probabilitate In cazul unui spatiu de esantionare discret, legea poate fi specificata dand probabilitati evenimentelor elementare, iar in cazul continuu dand probabilitati intervalelor 4Modelare si simulare Seminar 1 Aplicatia 2 Durata de prelucrare a unui apel telefonic variaza intre 0 si 1 secunde Considerand ca toate realizarile sunt echiprobabile, sa se determine: i) o lege de probabilitate potrivita pentru aplicatia analizate; ii) probabilitatea ca durata sa se indeparteze cu cel putin 0,3 secunde de centrul intervalului; Rezolvare: ≤ nu are sens, in i) Alocarea de probabilitati unei realizari ,(0 1)cc≤ este infinit nenumarabil acest caz, deoarece spatiul realizarilor S= si deci [{ }] 0Pc= Fie urmatoarea afirmatie: “Probabilitatea ca o realizare sa apartina unui subinterval al multimii S este egala cu lungimea intervalului” adica: [[ , ]]Pab b a=− pentru 01ab≤≤≤ Aceasta este o lege de probabilitate deoarece satisface axiomele ii) Evenimentul in discutie are urmatoarea forma analitica: A=∪ Intervalele fiind disjuncte putem aplica axioma 3 deci: [] 0,4PAP P=+ = Aplicatia 3 In urma masuratorilor timpilor de functionare ai unor componente electronice de acelasi tip s-a ajuns la concluzia ca: “proportia componentelor al caror timp ” Sa se determina legea functionare depaseste durata t scade exponential cu rata α de probabilitate care satisface acest rezultat Rezolvare: Spatiul de esantionare este: (0, )S=∞ Evenimentul continut in enunt {timpul de functionare depaseste durata t}A= =∞ se exprima analitic astfel: {;}(, )Axx t t=> Acestui eveniment i se aloca, conform afirmatiei, urmatoarea valoare: t−α [] [(, )]PAPt e=∞=, pentru 0t >0α> si Verificam ca aceasta alocare respecta axiomele 1, 3 t−α = ∀>α>⇒A i) 0[,] , 0, 0Pte t≤∞ 1 =⇒ ii) [] [(0, )] 1PSP A=∞ 2 5Modelare si simulare Seminar 1 ∞ iii) pentru evenimente de genul: (,Brs= si (, )Cs=, respectarea ] axiomei 3 impune: [( , )] [ , ] [( , )]Pr Prs Ps∞= + ∞ deci legea generala cautata este: ( ] ts−α−α [ , ] [( , )] [( , )]Prs P r P s e e=∞−∞=− ( ] Aplicatia 4 Fie un concentrator ce dispune de trei circuite de iesire Sa se stabileasca un model probabilistic al procesului de ocupare aleatorie a circuitele de iesire Rezolvare: Experimentul consta in ocuparea celor trei circuite de iesire Observatiile pot lua in considerare modul de ocupare (Varianta I) sau numarul circuitelor ocupate (Varianta II) Varianta I: Notand cu L starea de liber si cu 0 starea de ocupare, spatiul realizarilor este: { , 0, 0 , 00,0 ,0 0,00 ,000}S LLL LL L L L LL L L= I Putem presupune ca evenimentele sunt echiprobabile In acest caz fiecare eveniment elementar apare cu probabilitatea 1/8 Varianta II: Spatiul realizarilor este {0,1,2,3} Si in acest caz putem presupune ca realizarile sunt echiprobabile; Rezulta: fiecare eveniment elementar apare cu probabilitatea 1/4 OBS Desi corecte din punct de vedere teoretic, cele doua doua alocari de probabilitati conduc al valori diferite pentru probabilitatile acelorasi evenimente Astfel, in cazul evenimentului: {2 circuite din 3 sunt ocupate}A= ¾ pentru Var I obtinem: 3 [ ] PAPLLL== I 8 ¾ pentru Var II: 1 [] [{2}]PA P== II 4 Rezulta: Trebuie decis care este cea mai potrivita alocare in descrierea fenomenului 6Modelare si simulare Seminar 1 2 Probabilitati conditionate Probabilitatea conditionata ca evenimentul A sa se indeplineasca daca evenimentul B a avut loc este definita de relatia: []PAB∩ [/]PA B= []PB Aplicatia 5 La un nod de comutatie se receptioneaza mesaje ce sunt distribuite echiprobabil la un numar de N terminele identificate prin adresele 0,1, 2, , 1N− … Consideram urmatoarele evenimente: A = {mesajul receptionat se transmite unui terminal de adresa para} B = {mesajul receptionat se transmite unui terminal de adresa impara} C = {mesajul receptionat se transmite unui terminal de adresa M≤ unde } MN == ∑ 1kM=+ N Nk Nk− (1 ) 0,00259Cp p=−= k∑ 1kM=+ Pentru evitarea dificultatilor generate de termenii factoriali, in calculele computerizate se foloseste formula recursiva: ()nkp− (1)(pkpk)+= nn (1)(1 )kp+− (probabilitatea sa avem k succese in n incercari) Secventele de subexperimente independente in care se urmareste nr succeselor pentru mai multe evenimente reciproc exclusive sunt modelate cu ajutorul legii de probabilitate multinomiala Fie ,,,BB B… evenimente ce realizeaza o partitie a spatiului de esantionare S 12M Probabilitatile de aparitie a lor in cadrul unui subexperiment, []PBp= , satisfac jj + =… relatia: 1pp p++ 12M Fie k numarul de succese ale evenimentelor ,1,,BjM= …, obtinut in n jj ce specifica numarul succeselor repetitii independente Vectorul (, , , )kk k… 12M pentru fiecare eveniment satisface legea de probabilitate multinomiala: 10Modelare si simulare Seminar 1 !kk kn 12M [( , , , )]Pkk kp p p=⋅ ……⋅⋅ 1212MM !! !kk k⋅⋅ ⋅… 11M kk k++nunde: + =… 12M Aplicatia 8 Un calculator multitasking, multiuser executa 3 categorii de programe caracterizate prin duratele 10msec, 100msec si 200msec Aceste programe sunt cerute de utilizatori cu probabilitatile: 0,2; 0,3 si 0,5 Sa se determine probabilitatea executiei simultane a cate 3 programe din fiecare categorie Rezolvare: Multimea programelor de orice durata S se partitioneaza in: multimea progde 10msec cu [ ] 0, 2;BP== B 11 B multimea progde 100msec cu [ ] 0,3;BP== 22 multimea progde 200msec cu [ ] 0,5;BP== B 33 Evenimentul A = {se executa simultan cate trei programe din fiecare categorie} are expresia analitica: A=(3,3,3) Aplicand legea de repartitie multinomiala avem: 3339! [(3,3,3)] (0, 2) (0,3) (0,5) 0,04536P =⋅⋅= … 3!3!3! Secventele de subexperimente Bernoulli in care se urmareste numarul de repetari pana la aparitia primului succes sunt modelate cu ajutorul legii de probabilitate geometrica 1CC Cm− [& && ](1) , 1,2,Pm PA A A A p p m==−⋅……= () 12 1mm− unde m este numarul de repetari pana la aparitia succesului CC C ¾ ,,AAA… este secventa de insuccese obtinute in primele m-1 incercari; 112m− ¾ A este succesul evenimentului in incercarea m; m ¾ p este probabilitatea aparitiei succesului intr-o incercare; Probabilitatea ca numarul necesar de incercari m, pentru obtinerea succesului sa fie k − mai mare decat un numar fixat k este: [](1)Pmk p>= 11